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Abstract 

We consider CP- violating effects in r — > Ktitiv T) assuming that a charged Higgs 
boson provides a new amplitude that can interfere with the usual Standard Model 
amplitude. We consider four CP-odd observables - the regular rate asymmetry, 
two modified rate asymmetries and a triple-product asymmetry. The regular rate 
asymmetry is expected to be small because it requires the interference of the new 
physics amplitude with the standard model amplitude containing the hadronic scalar 
form factor. The other CP asymmetries may be more promising in terms of their new 
physics reach. Numerical estimates indicate that the maximum obtainable values 
for the modified and triple-product asymmetries are on the order of a percent. 
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1 Introduction 



In the standard model (SM) of particle physics, CP violation is due to a complex 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. But is this the only 
source of CP violation? In order to answer this question, it is important to look for 
CP-violating effects in as many systems as possible. 

One such system is r decays. In the SM, CP violation in the r system is essen- 
tially zero pQ; we consider, instead, a search for physics beyond the SM. In Ref. [2], 
we examined CP violation in strangeness-conserving r decays. It is only natural next 
to turn to those processes with AS = 1. The simplest such decay is r — ► Ki\v T . 
However, CP violation in this process has been analyzed in detail in Ref. [3], and we 
have nothing to add here. The next decay is r — > K-k-kv t . This has been examined 
theoretically in the past in Refs. [H In this paper we update these analyses [6J. 

One has to assume the presence of new physics in order to get non-zero CP- 
violating effects when comparing t~ — > K~7r~7T + u T to its CP-conjugate decay. In 
Ref. [I], the left-right (LR) model is assumed when the authors consider r — > Kiriri/ T . 
However, as shown in Ref. j2], if there is no LR mixing, CP violation is proportional 
to the mass of the neutrino, and is negligible. Thus, in the LR model, CP violation 
in t decays is proportional to Wl-Wr mixing. However, we know this is small [7]. 
We therefore conclude that sizeable CP violation in the r system cannot arise in the 
LR model. 

For this reason, in this paper, we assume that the r decay includes the exchange 
of a new-physics (NP) charged Higgs. Note that many NP models have two Higgs 
doublets, so that a charged Higgs is present. However, if the Higgs doublets give mass 
to the fermions, the coupling of the charged Higgs boson is generally proportional 
to the masses of the first- and second-generation quarks. Since these are small, CP 
violation in the r system will also be small. To avoid this, if CP violation is to be 
observed in r decays, the charged-Higgs coupling must be large. In other words, 
t — ► KttttUj- probes non- "standard" NP CP violation. 

It is worth noting at this point that CLEO has searched for CP violation in 
r — > Knu T [S] and has set a bound on a coupling constant related to the scalar 
coupling of a charged Higgs (or other scalar boson) to the up and strange quarks. 
The experimental investigation suggested in this work would be complementary to 
that carried out in Ref. [8] in that it would probe the pseudoscalar coupling of the 
Higgs to the up and strange quarks. In the notation introduced below, the CLEO 
experiment probed T]s, while a CP analysis of r — > K-R-nv T would probe r]p [see 
Eq. Q below]. 

In the presence of one NP contribution, the amplitude for the decay r — > Ktt7tv t 
can be written 

A = A x + A 2 e^e lS , (1) 

where <fr and 5 are the relative weak (CP-odd) and strong (CP-even) phases, respec- 
tively. The amplitude for the antiprocess, A, is given by the same expression, but 
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with cj) — > —(f). 

In general, CP violation is obtained by comparing \A\ 2 to \A\ 2 . There are three 
types of signals: 

1. The full rate for a particular process involves Dspms l-4| 2 ; integrated over the 
final-state momenta in the usual way. The rate asymmetry is given by the rate 
difference of the process and antiprocess. 

2. The rate asymmetry can be altered in two ways. First, if some spins are mea- 
sured, one does not sum over them. Alternatively, one can integrate asym- 
metrically in order to isolate certain terms in the differential width. In either 
case the process-antiprocess difference leads to a modified rate asymmetry. 

3. One can also construct CP asymmetries based on the quantity V\ ■ (v 2 X v 3 ), 
where each Vi is a spin or momentum. This is a triple product (TP), and its 
value can be different for process and antiprocess, signaling CP violation. 

The rate asymmetry or modified rate asymmetry is proportional to 

sin S sin (2) 

(integrated over phase space). Thus, this category of CP violation requires that 
the two decay amplitudes have a non-zero relative weak and strong phase. The TP 
asymmetry is proportional to 

cos 5 sin , (3) 

so that one does not require a strong-phase difference to get a TP asymmetry. In 
this paper we consider all three types of CP violation in r — > Kintv T . (Ref. [I] 
considers only TP's.) 

The remainder of this paper is organized as follows. In Sec. [2j we write down 
the expression for the differential width for r~ — > K~7T~7T + b> T in terms of various 
form factors, and including the NP contribution. We perform weighted integrations 
of the differential width over phase space to isolate certain cross-terms. In Sec. [3j 
we consider four CP asymmetries: the regular rate asymmetry, two modified rate 
asymmetries, and a triple-product asymmetry. The modified rate asymmetries and 
triple-product asymmetry are constructed using the weighted differential widths 
from Sec. [2] Section [4] contains a numerical analysis of the modified rate asymmetries 
and the triple-product asymmetry. We conclude with a few closing remarks in Sec. [5j 

2 Differential width for r~ — > K~n~7T + u T 

We start by determining an expression for the differential width for t~ — ► K~tt~tt + v t , 
including possible NP effects due to a new charged Higgs boson H . 
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2.1 General expression for the differential width 

Let us begin by considering the SM contribution to t~ — > K~7r~7T + v T . Within the 
SM, the relevant effective Hamiltonian is given by 

ftfM = ^ sin Q c V T% {1 - 75 )r s^(l - lb )u + h.c, (4) 

where 8 C is the Cabibbo angle. The hadronic matrix element for the decay may be 
conveniently parameterized in terms of four form factors as follows [5], 



J" = (K-(p 1 )7r-(p 2 )7r+(p 3 )|s7 /i (l-7 5 )w|0) 

F 1 (s 1 ,s 2 ,Q 2 )(p 1 - p 3 ) u + F 2 (si, s 2 ,Q 2 )(p 2 -P3 



+iF 3 (s u 8 3 , Q 2 y up °p lv p 2p p 3(T + F 4 (si, s 2 , g 2 )Q" , (5) 

where = {pi+p 2 +p 3 Y, = g^-Q»Q v /Q 2 , s 1 = (p 2 +p 3 ) 2 and s 2 = (pi+p 3 ) 2 ; 
also, we adopt the convention eoi23 = +1 as in Ref. [9j Q The form factors F1-F4 have 
been considered, for example, in Ref. [10J. As noted there, F% can arise due to the 
decay chain r — > K\V T , with K\ — > K*ir — > Ktttc, while F 2 comes from r — > KiV T , 
with Kx-^Kp-^ Kirn. It is now known that both the ^(1270) and the ^(1400) 



contribute (see Sec. 4.1 for further details). F 3 is the anomalous Wess-Zumino term. 
It can be estimated by considering the decay chain r — ► K*i/ T , with the intermediate 
K* going to pK or K*n [10J. The scalar term, F 4 , is generally assumed to be 
negligible for this decay, since there is no suitable pseudoscalar resonance through 
which the decay can proceed. The authors of Ref. [H] performed a calculation of 
F4 within the context of Chiral Perturbation Theory and found that F4 is non-zero 
if one includes chiral-symmetry-breaking mass terms for the quarks. The resulting 
expression for F4 was found to contain both a non-resonant term (proportional to 
m 2 + m 2 K ) and a resonant term. A numerical study indicated that the SM scalar 
contribution to the width was quite small [11]. We will consider the form factors 



further in Sec. 4.1 At this point we simply note that Fi and F 2 give the dominant 
contributions to the rate for r — > K-n-nv T [12] , while numerical estimates indicate 
that the Wess-Zumino term (F 3 ) gives a subdominant contribution. In fact, in 
their experimental analysis, CLEO discards the Wess-Zumino term altogether and 
considers only the contributions due to F\ and F 2 [T2] . 

Starting from Eq. (Jsl), the amplitude squared for t~ — > K^-n^-n + v T within the 
context of the SM is given by 

I^smI 2 = ^sm 2 e c L^\ (6) 



lr The authors of Ref. [9] adopt the convention £0123 = +1, but don't state the precise functional 
form for F 3 . Subsequent authors state F 3 , but the sign of eoi23 is not obvious. We make a particular 
choice for the sign of F3 below; changing this sign would change the sign of the related asymmetry. 
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where = M M (M u y and = J" (r)\ with M M = u Ut ^(1 - 7 5 )m t . 

Effects due to a charged Higgs modify the effective Hamiltonian relevant for 
t~ — * K~7r~7t + u T , adding the following termsj^] 

G 

H-ei = sin0 c [ri s v T {l + 75 )r su + r) P v T {l + 75 )r s 75 w] + h.c. (7) 

The total effective Hamiltonian is then H e s = "Hfjj 1 + 'H^z- In writing down Eq. ^ 
we have neglected terms that would involve a right-handed projection of the neutrino 
field. The interference of such terms with the SM amplitude would be suppressed 
by m Ui (assuming that the neutrino spin states are summed over). 

The NP effects can be incorporated into the amplitude in a straightforward man- 
ner. We first define a new current J M , which is obtained from J M by the replacement 

F 4 P 4 = F4 + JjL VPj ( 8 ) 
m T 

where the pseudoscalar form factor has been defined as follows 

(K-( Pl )n-(p 2 )n+(p 3 )\s^ U \0) = f H . (9) 

Defining H^ u = J M (j v ^j , we then find the following expression for the square of 
the matrix element, 

\A\ 2 = ^sin^V^. (10) 

Note that we have used the r~ equation of motion in order to arrive at our definition 
of F4. We have also neglected the mass of the neutrino. 

The decays of r leptons to final states containing two and three pseudoscalar 
mesons have been thoroughly analyzed in Ref. [9]. The notation described there 
has been adopted widely in the field and is quite standard. First, let us define 
several useful angles. Our definitions are identical to those in Ref. [9]. We review 
the various definitions here for convenience (more details may be found in Ref. [9]). 
The angle 6 is defined in the r rest frame; in that frame it is the angle between 
the direction of the hadrons ("Q") and the direction of the tau in the laboratory 
frame. All other angles are defined in the hadronic rest frame (i.e., the frame in 
which Q = pi + f>2 + f>3 = 0). In the hadronic rest frame we define two different 
coordinate systems, S and S'. These two coordinate systems are related by an Euler 
rotation using the Euler angles a, [3 and 7, as indicated in Fig. [Tj The z' axis in S' 
is chosen as the direction of the laboratory in the hadronic rest frame (n^). The x' 
axis is chosen such that the r direction (n T ) is in the x'-z' plane, making an angle 

2 These expressions are similar to those in Ref. [4J, although our notation differs slightly from 
that found there. 
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A 




Figure 1: Definitions of the angles a, /3, 7 and ip. The Euler rotations corresponding 
to a, (3 and 7 are about the z', y\ = 1/2 and z axes, respectively. This figure is very 
similar to a figure found in Ref. [9]. 



ip with respect to the z' axis (see Fig. Q. The z axis in S is perpendicular to the 
plane defined by the momenta of the hadrons: z = h± = pi x P2/ \Pi x p* 2 1 - The x 
axis is taken to be the direction of p 3 ; i.e., x = p 3 / \p 3 \. The three Euler angles are 
defined as follows: a is the angle between the (n^, h T ) plane and the (n^, h±) plane, 
j3 is the angle between hi and h± and 7 is the angle between the (fiL,h±) plane and 
the x) plane. 

Having defined the various angles, we may write the differential width for t~ — > 
K~7t~tt + v t as follows [9], 



G^sin 2 ^ m 2 T -Q 2 dQ 2 da d-y d cos (3 d cos - 

al = s — 7t ctSi aS9 L,, v ri^ , (11, 

256(2vr) 5 m T m 2 Q 2 2tt 2tt 2 2 ' V y 



where Q 2 , S\ and s 2 were defined below Eq. ([5]). 

The coordinate system 5* is convenient for expressing the momenta of the three 
pseudoscalar mesons and for computing the various components of the tensor H^ u . 
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In this coordinate system we have [9], 



Pi 
/' 

Pi 



where 



Ei = (Q 2 - Si + m 2 ) / (2^) 



Pi 
p\ 



Pi 



v 
Pi 

u 

pi 



El 



m 



3 • 



(2E 1 E 3 -s 2 + m\ + ml) / {2p 3 ) 
(2E 2 E 3 - Sl +ml + mfj / (2p x 3 ) 



El - {pi) 



m 



i i 



-\JEl - (pg) -m 2 2 = -p\. 



(12) 
(13) 
(14) 



(15) 

(16) 
(17) 
(18) 
(19) 
(20) 



S3 is defined analogously to s\ and s 2 [i.e., S3 = (p± + P2) 2 ] and may be expressed in 
terms of S\, s 2 and Q 2 . Note that the angle between p\ and p 3 is fixed for a given 
choice of S\, s 2 and Q 2 . 

The above definitions for the^various hadron momentum vectors allow us to 
determine simple expressions for H^ u in S. We will not write out all 16 elements 
of the tensor. Rather, we define new quantities Bi (i — 1, . . . , 4) that are related to 
the components of J M as follows, 



B 1 
B 2 

B 3 
B A 



J^[F 1 (p 1 -p 3 ) x + F 2 (p 2 -p 3 ) x ] 
J 2 = (F 1 -F 2 )p», 

-u 3 = f 3 Jq^pIp x , 



m 1 f H 

Fa H Vp 

m T 



(21) 
(22) 
(23) 

(24) 



Then the components of consist of various combinations BiB* (in some cases 
multiplied by ±i). Inserting these expressions into Eq. ( JTT[) and integrating over a 
we find an expression very similar to that given in Ref. [4JPJ 



dY 



G 2 F sm 2 6 c (m 2 T -Q 2 Y 



dQ 2 dsi ds 2 d'y d cos (3 d cos 6 512 (27r) m\Q 2 



\k 1 + K 2 + ^K 1 (3cos 2 /3-l)/2 



Bi\ 



I -Bo 



3 Due to some ambiguities, it is difficult to tell if the expressions agree exactly. 
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\k x +K 2 -\k x (3cos 2 /3-l)/2 



\B 3 \ 2 + K 2 \BA 2 



K x sin 2 /3cos 27 ( | Si | 2 - |P 2 | 2 ) + K x sin 2 /? sin 2 7 Re(B x B. 



+ 



+ 2K 3 sin /3 sin 7 Re(S XJ B*) + 2i? 2 sin /3 cos 7 Re (Pi B* 4 ) 
+ 2F3 sin /3 cos 7 Re(P 2 P*) - 2Z 2 sin (3 sin 7 Re(P 2 P^ 
+ 2K 3 cos /3 Im(5 x ^) + K x sin 2/3 cos 7 Im^Pg*) 

- sin 2(3 sin 7 Im^K:) + 2F 2 cos (3 \m(B 3 Bl) 



(25) 



Note that, of the four parameters Bi defined in Eqs. (21)-(24), only P4 contains a 



non-SM weak phase. Thus, the only terms in Eq. (25) that can lead to non-zero 
CP asymmetries are those containing one or more powers of P4. The parameters 
Ki and Ki in the above expression are defined as follows [9], 



K x 


= 1 -Pcosfl- (m 2 /Q 2 ) (1 + Pcosfl), 


(26) 


K 2 


= (m 2 /Q 2 ) (1+P cos 9), 


(27) 


K 3 


= l-Pcos9, 


(28) 


K x 


= K 1 (3cos 2 ^-l)/2-(3/2)ir 4 sin2^, 


(29) 


K 2 


= K 2 cos ip + K^smip , 


(30) 


K 3 


= K 3 cos if) — K4 sin if) , 


(31) 


K 4 


= Jm 2 /Q 2 P sin 9, 


(32) 



where the parameter P denotes the polarization of the r~, s 2 = — P 2 . In the nu- 
merical work in Ref. jl], a value of P was used that was relevant for LEP. In our 
numerical work we will take P = 0, which is appropriate for lower-energy experi- 
ments [9]. Note that if> is a function of cos# and Q 2 . If the r's are pair-produced at 
a symmetric collider, 



cos 9 



cosV' 



(2xm 2 —m 2 — Q z 



{m 2 - Q 2 ) sjl - Am 2 /s 

x (m 2 + Q 2 ) - 2Q 2 
(m? - Q 2 ) Jx 2 - AQ 2 /s 



(33) 
(34) 



where x = 2E\ l j\fs and s = 4P 2 eam , with Eh being the hadron energy in the lab [9] 
(see also Ref. [13]). Thus, given s (we take s = (10.58 GeV) 2 below), Q 2 and cos#, 
one can solve for x and substitute this expression into the expression for cosip. 
Finally, note that if the direction of the r~ could be determined, it would not be 
necessary to integrate over a. In this case it might be possible to extract other 
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useful information for the construction of CP asymmetries. We do not consider this 
possibility in this work. 



The differential width in Eq. (25) may now be integrated to compute the partial 
width for r~ — > K^-k^-k + v t . Comparison with the analogous quantity for the r + 
decay yields the regular rate asymmetry. One can also integrate over the angular 
variables in an asymmetric manner in such a way that certain cross terms are selected 



from Eq. (25). These "weighted differential widths" can then be compared to the 
analogous expressions for the r + decay to yield CP-odd quantities. We consider two 
types of asymmetries formed in this manner - modified rate asymmetries [whose 
dependence on the strong and weak phases is given in Eq. ^2j)] and a triple product 
asymmetry [see Eq. (131)]. 



2.2 Weighted differential widths 

The authors of Ref. pi] derived an expression for the differential width that is very 



similar to Eq. (25). Since they assumed that /# = for these decays, they only 
considered LR effects. We consider the complementary point of view. Assuming that 
fn could be non-zero and noting that there are strong constraints on LR mixing, we 
consider only effects due to the exchange of a charged scalar. The analysis in Ref. [I] 
focused exclusively on triple products in the differential width. In our notation, these 
TP's correspond to the cross-terms containing the factors lm(BiB*). Recall that 
CP asymmetries formed from triple products do not require the presence of strong 
phases [see Eq. (|3j)]. 

In this work, we reconsider CP asymmetries formed from triple products and also 
consider CP asymmetries that can be formed from T-eveE0 cross-terms in the differ- 
ential width. Both types of terms may be isolated by employing suitable weighting 
functions when performing the angular integrations. 

We begin by defining various regions in terms of 7 and (3, as in Ref. [I], 

I : < 7 < 7r/2, II : vr/2 < 7 < vr, 

III : vr < 7 < 3vr/2, IV : 3vr/2 < 7 < 2tt; , . 

(35) 

A : < (3 < tt/2, B : vr/2 < (3 < vr; 



As noted above, in order for a particular term in the differential width [Eq. ( [25] )] 
to contribute to a non-zero CP asymmetry, it must contain one or more powers 
of -B4. This is because B4 contains the possible CP- violating phase coming from 



NP. Inspection of Eq. (25) leads one to the conclusion that there are four terms of 
interest. One is proportional to I-S4I . As we shall see below, this term arises in the 
regular rate asymmetry. The remaining three terms are proportional to the angular 
functions sin (3 cos 7, sin/? sin 7 and cos/3. These three terms can be isolated by 

4 "T-even" here refers to the naive time-reversal operation. 



S 



Table 1: Angular weighting factors. The regions I-IV, A and B are defined in 
Eq. (35) in the text. The second column gives the angular functions of interest, 
/i(7, 0). The third column gives the weighting function 5^(7, 0) that can be used to 
isolate f%( r y,0)- The functions gi{ r y,0) are simply ±1 depending on which region 7 
and fall in. 



i 


Ml, (3) 


9i(l,0) 




1 


sin sin 7 


I + II - III - IV; A 


+ B 


2 


sin cos 7 


I - II - III + IV; A 


+ B 


3 


cos 


I + II + III + IV; A 


-B 



using appropriate weighting functions, as indicated in Table [TJ Thus, for example, 
to isolate the term in Eq. (25) proportional to sin sin 7, the differential width is 
multiplied by ^1(7, 0) (which is +1 in regions IA, IIA, IB and IIB and —1 in the 
other regions) and the angular integration is carried out. This eliminates all other 
terms since the weighting functions are such thatj^] 



fi{l,0)9j{l,0) sin 0dqd0 



2tc 5ij 



1,2,3). 



(36) 



Using the weighting functions g%{pi,0) in Table [TJ we define weighted differential 
widths as follows, 



dQ 2 dsi ds$ 



dT 



dQ 2 dsi ds 2 c?7 d cos d cos 9 
The results for the three weighting functions are as follows, 
dT x 



<7i(7, 0) sin d0 dj dcos9 . (37) 



dQ 2 ds\ ds2 

dY 2 
dQ 2 dsi d$2 

dV 3 
dQ 2 ds\ ds 2 



A(Q 2 ) 
A{Q 2 ) 
A(Q 2 ) 



where 



A(Q' 



{K 3 )Re(B 1 B* 3 ) - (K 2 }Re(B 2 Bl) 
(K 3 )Re(B 2 B;) + (K^ReiB.Bl) 
(K 3 )Im(B 1 B*) + (K 2 )lm(B 3 Bl) 

G 2 F sin 2 9 c K-Q 2 ) 2 
128 (2tt) 5 m\Q 2 ' 



(38) 
(39) 
(40) 

(41) 



5 The weighting functions are also orthogonal to 1, (3 cos 2 /? — l) /2, sin 2 /3 cos 27, etc., so that 
only the intended cross-terms are isolated. Also note that, experimentally, a more statistically 
significant weighting procedure might be to weight the differential width by the various functional 
forms fi themselves. See Ref. [T3] and also the moment analysis discussion in Ref. [9]. 
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(Ki) = 



Kisinede. 



(42) 



The three weighted differential widths defined in Eqs. (38)-(40) can now be 
compared to their CP-conjugates in order to construct CP asymmetries. Recalling 
that the Higgs contribution resides in i? 4 [see Eq. (24)] and noting that each of 
the three expressions above contains a term linear in B4, we see that each of the 
resulting CP asymmetries has the possibility of being non-zero. 

In the following sections we construct the CP asymmetries and then study them 
numerically to see if they might provide useful probes of non-SM CP violation. 



3 CP-odd observables 

Before analyzing the various CP asymmetries, let us consider the coefficients Bi 



defined in Eqs. (21)-(24) a bit more carefully. As noted above, the sole non-SM 
weak phase resides in _B 4 . The form factors Fj and /# are potential sources of 
strong phases. We may thus parameterize the four coefficients as follows, 

(43) 
(44) 
(45) 
(46) 



By = 


filled, 


B 2 = 


B 2 \e i5 \ 


B 3 = 


B 3 \e lS \ 


B 4 = 


bP 


e iSi 



where 



B 



(i) 



! Q 2 K 



■1 ■ 



B 



+ 



(2) 



B 



(2) 



i5 H +i<j>i 



'Q 2 —Vp, 



(47) 



and where 5i and <$h represent strong and weak phases, respectively. An explicit 
expression for the weak phase 4>h is as follows, 

4> H = arg (r] P ) . (48) 

This phase could in principle be of order unity. 

As was the case in Ref. [2], we can consider three types of CP asymmetries. 
The first is the regular rate asymmetry. This asymmetry is likely to be small in 
t — > K7i7iu T and is therefore unlikely to be measureable in the near future. The 
second and third types of asymmetries are the modified rate asymmetry and the 
triple-product asymmetry. We consider two different modified rate asymmetries, 
and one triple-product asymmetry. The triple-product asymmetry is similar to one 
considered for the decay r — > KttKv t in Ref. [3]. The modified rate asymmetries, 
to our knowledge, are new relative to this decay mode. Both types of asymmetries 
are constructed by first performing an asymmetrical integration over the kinemat- 
ical angles j3 and 7, as noted in Eq. (37) and Table [TJ Since the procedures for 
extracting these two types of asymmetries are similar, we consider them together in 
the following. 
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3.1 Rate asymmetry 

Let us first consider the regular rate asymmetry. In this case the angular integrations 
are performed symmetrically [#(7, 0) = 1] and the width for the process is compared 
to that for the anti-process. The differential width for the r~ decay in this case is 
given by 



dT 



dQ 2 dsi ds^ 



A{Q 2 ) 




(49) 



The width for the r + process will have the same strong phases, but the weak phases 



will have their signs reversed. It is immediately evident from Eq. (49) and Eqs. (43)- 



(45) that the coefficients B\, B2 and B% will not give any contribution to the rate 



asymmetry, since they do not contain weak phases. Thus, the rate asymmetry is 
proportional to^] 



\Ba 



B, 



4 



B 



(i) 



B 



(2) 



sin (64 - 6 H ) sin ((j) H ) 



(50) 



This expression is proportional to \F i fnrip\. The SM scalar form factor F 4 is gen- 
erally thought to be small. If the NP factor r]p is also small, then the regular rate 
asymmetry is doubly suppressed. Given the expected smallness of the regular rate 
asymmetry, we do not consider it further here. As we shall see, however, other 
CP asymmetries can be constructed that depend on FifnT/p, with i — 1, 2, 3. Such 
asymmetries may be more promising in terms of their NP reach. 



3.2 Modified and triple-product CP asymmetries 

We define CP asymmetries corresponding to the weighted differential widths [Eqs. 
(PI-pOll as follows, 



A 





CP 



dT; 



dT,, 



T + T 



dQ 2 ds\ ds2 dQ 2 ds\ ds2 



dQ 2 dsi ds2 



(51) 



The quantities with the bars correspond to the decay r + — > K + 7i + ti v t and are 
obtained from those without the bars by changing the signs of all weak phases 



while leaving strong phases unchanged 



A^p and A^'p descend from the terms 



I (2) 



6 This expression is part of an integral over phase space. Note that one or both of the strong 
phases could depend on Q 2 , si and S2- 

7 Note that we subtract the width for the anti-process from that for the process, both for the 
modified rate asymmetries and for the triple-product asymmetry. The authors of Ref. @] consider 
only triple-product asymmetries. Their expressions for the anti-process contain an extra over-all 
sign; thus they add the widths for the process and anti-process to obtain CP asymmetries. This 
is a notational difference. Both approaches lead (correctly) to a TP CP asymmetry that is of the 
form of Eq. pi). 
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containing Re(B2Bl) and Re(BiBl) in Eqs. (38) and (39), respectively. These are 

(3) 

both modified rate asymmetries. The third asymmetry, A CP , descends from the 
term containing Im^B^Bl) in Eq. (40). This a triple-product asymmetry. T and T 

K~7T~7T + h> T and r + 



in Eq. (51 ) represent the partial widths for r 



K 



IT ' IT V 



respectively. In our numerical work below we make the approximation that T 
so that T + T ~ 2r. The experimental value for T is used. 



r. 



3.2.1 Modified rate asymmetries (i = 1,2) 

The modified rate asymmetries, A^ P and A^ P , require a strong phase in order 
to be non-zero. These asymmetries are analogous to the "polarization-dependent 
asymmetry" defined in Ref. [2]. In order to obtain numerical estimates for these 
asymmetries, let us make the following simplifying assumptions. First of all, we will 
assume that fn has no Q 2 , Si or s 2 dependence. We will also assume that fn has 
no strong phase associated with it (it will be taken to be real and positive). Under 
these assumptions, these two asymmetries are given by 



.4 



(i) 

CP 



r + r 



A(Q 2 



cosipVi Im(-Fi — F 2 ) dQ 2 dsids 2 dcos9 



x f H Im(r) P ) 



(52) 



.4 



(2) 
CP 



r + r 



A{Q 2 



x f H lm(rj P ) 



cos ip Im[Fi (p 1 - p 3 ) x + F 2 (p 2 - Ps) x ] dQ 2 ds 1 ds 2 dcos 9 



(53) 



in which we have taken the r's to be unpolarized (P = 0). Recall that ip depends 



on 9 through Eqs. (33) and (34). 



As noted above, A^ P and A^ P both have the generic form sin 
Im(r) P ) oc sm<j) H and lm(F 1 - F 2 ) and Im[Fi (p% - p 3 ) x + F 2 (p 2 - p 3 ) x 



sin 5, since 
are both 

proportional to sin 5, with 5 being a strong phase. In Sec. |4.2| we will examine the 
sensitivity of these asymmetries in a particular model for the form factors. 



3.2.2 Triple-product asymmetry (i = 3) 

The third CF asymmetry, A CP , is a triple-product asymmetry and is similar in some 
respects to the asymmetries constructed for r — > Kttttv t in Ref. [1] . Recall, however, 
that in that case the authors assumed that the NF effects were due to a new right- 
handed gauge boson. To obtain a numerical estimate for A CP , we make the same 
simplifying assumptions as above; i.e., we assume that fu is real and positive (no 
strong phase) and that it has no Q 2 , Si or s 2 dependence. Under these assumptions, 



.4 



(3) 
CP 



r + r 



A(Q 2 ) cos ip p\ pi Re(F 3 ) dQ 2 ds 1 ds 2 dcos 9 



12 



x f H lm(rj P ) 



(54) 



Like the modified rate asymmetries considered above, the triple-product asymmetry 
A^p is proportional to Im(?7p). In contrast to A^ P and Aq P , however, this asym- 
metry does not require a strong phase, since Re(_F 3 ) oc cos S (where 5 represents a 
strong phase). Having said this, there is a potential drawback with A^ P in that it 
depends on the sub-dominant Wess-Zumino form factor F 3 , whereas A^ P and A^ P 
depend on combinations of the dominant form factors F\ and F 2 . In the next section 
we perform a numerical study to examine these various factors quantitatively. 



4 Numerical Results 

The modified and triple-product asymmetries defined above all have the form 

4p = a^pfnlmM , 2 = 1-3, (55) 

where the constants determined by integrating over cos^, s\, s 2 and Q . 

In this section we assume particular functional forms for the form factors and use 
these to estimate the cl§p- ^ turns out that there are significant cancellations that 
occur as one performs the integrations over phase space. To help illustrate this 
cancellation, we define four differential quantities as follows, 



da 



(*) 

CP 



dX 



(56) 



with X given by M K7T7T = yQ 1 , M W7T = y/si, M Kw = and cos 8. Given the 
cancellations that occur upon integration, experimentalists may wish to study dif- 
ferential CP asymmetries in addition to, or in place of, the integrated asymmetries. 



4.1 Model for the form factors 

There have been several models for the form factors F1-F3 over the past number 
of years. One model, which simply took the intermediate K\ to be the ^(1400), 
may be found in Ref. [HI] (see also Ref. [IS])- A subsequent analysis by Finkemeier 
and Mirkes [18J took into account both the i^i(1400) and the ^(1270) resonances 
and also incorporated other K* resonances {K*' and K*"; see also Refs. P22 120]). 
Finally, an experimental analysis of the form factors was performed by the CLEO 
collaboration in Ref. [12] . 

The various models that have been proposed make different assumptions regard- 
ing the anomalous Wess-Zumino term, F3. The authors of Ref. pj0] found that the 
F 3 term contributes approximately 1% to the overall width for t~ — ► K~7t~7t + i> t . 
The parameterization in Ref. JTS] led to an anomalous contribution of order 10%. 
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The CLEO collaboration noted that the contribution would be of order 5.5% based 
on a particular model (found in Ref. [15]). Since the contribution was expected 
to be small, they set F 3 to zero in their analysis and focused on determining the 
resonance structures of F\ and F 2 . The uncertainty resulting from the neglect of F 3 
was incorporated into their systematic error |12j . 

We model our numerical work after the CLEO analysis, with the main exception 
being that we allow F 3 to be non-zero. Guided by Ref. [12] for Fi and F 2 and by 
Ref. [TO] for F 3 , we write the form factors in terms of various Breit-Wigner functions 
as follows, 



F 1 (s l ,s 2 ,Q 2 ) 
F 2 {s u s 2 ,Q 2 ) 

F 3 (si,s 2 ,Q 2 ) 



2N 



3F n 


C 


N 




V3F 




N 3 





1270 



(Q 2 ) + D- BW imo {Q 2 ) BW K *(s 2 ), (57) 



2y/2n 2 F* 



A ■ BW l2m (Q 2 ) + B ■ BWuooiQ 2 )] T^{ Sl ) , (58 



BW K *{Q 2 ) x 



TM{s{) + clBW k .{s 2 ) 
1 + a 



(59) 



with a = —0.2 and F n = 93.3 MeV. Also, we set F4 to zero and only take fn 
into account when computing the numerators of the asymmetry expressions. The 
constants N, N 3 and A-D will be discussed further below. The normalized Breit- 
Wigner propagators for the lfi(1270) and the Ki(1400) are assumed to be given 
by 



bw k M 



-m 2 Ki +im Kl T Kl 
Q 2 -m 2 Ki +im Kl T Kl ' 



(60) 



with rriKi and being the mass and width for the appropriate K% state. As noted 
in the CLEO analysis, a fit to the r — > Kt^txv t data indicates that the effective 
i^i(1270) and Z£i(1400) widths are larger in this decay than the respective values 
reported by the Particle Data Group (see also Refs. [16, 17]). Following CLEO, we 
take the following values for our numerical analysis [12] : 

mi 270 (i4oo) = 1-254 (1.463) GeV, r 1270 ( 14 oo) = 0.26 (0.30) GeV. (61) 

The Breit-Wigner propagators for the K* and p are taken to have energy-dependent 
widths (see, for example, Refs. [T2~l ITU] ). 



BW R {s) 



-m 



R 



with 



V R (s] 



s -m 2 R + iy/sT R (s) 



s \Pr) 



(62) 



(63) 
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where 



P = 7Tr\J[ s ~ ( m i +m 2 ) 2 }[s - (m x - m 2 ) 2 ] , (64) 

l\J s 

Pr = ^— \J[m 2 R - (mi + m 2 ) 2 ][m R - (m x - m 2 ) 2 ] . (65) 
When using the above expressions it is assumed that the resonance R decays to two 



particles with masses mi and m-i- [Equation (63) also assumes that \fs > mi + m 2 
- otherwise r#(s) should be set to zero. This condition is satisfied in all regions of 
phase space for the decay chains that we consider.] For the K*, a single resonance 
(with an energy-dependent width) is assumed; we take m K * = 0.892 GeV and 
TV* = 0.050 GeVj^The expression for the p incorporates two different resonances 
(the p and the p'), 

r« M = WM + WM) , (66) 

with (3 = -0.145, m p = 0.773 GeV, m p , = 1.370 GeV, T p = 0.145 GeV and 
= 0.510 GeV HBl El. 



Let us now consider the constants N, N 3 and A-D in Eqs. (57)- (59). The CLEO 
collaboration effectively set iv^ = in their analysis and then determined A-D [12J. 
The overall normalization N was not stated. We choose values that are similar to 
those reported in Table I of Ref . [32] , 




6917 0.195 



61636 

0.266 . (67) 

An apparent typo in Eq. (2) of Ref. [12] renders the relative signs of the constants a 
bit uncertain. The signs we have chosen for A-D are consistent with the signs used 
in Ref. [H] . Our parameter choice gives results for the differential width (see Fig. [2] 
below) that are visually similar to the results obtained in Ref. [12], although the 
agreement between our numerical results and those of CLEO is not perfect. Since 
we wish, in part, to study effects due to the inclusion of the F 3 term, we retain 
a non-zero value for AT 3 . As was noted above, there have been various estimates 
regarding the F3 contribution to the width, with estimates varying from 1% to 10% 
in papers that we have noted. For the purpose of our numerical study, we fix N and 



8 Note that the intermediate K* represents a K*° in the expression for the form factor Fi, while 
both K*° and K*~ appear in F 3 . For simplicity we use the same mass and width for both the 
charged and neutral versions of this particle. 
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Figure 2: Plots of the differential widths dT/dM, including the contributions from 
the various decay chains. The p, K*, i^i(1270) and i^i(1400) curves include contri- 
butions from only the -F\ and F 2 terms. The "W-Z" curves represent the contribution 
from the anomalous Wess-Zumino form factor, F 3 . 



N3 such that the F3 term contributes 5% to the r — > K7nri> T width, with the F\ and 
F2 terms contributing the remaining 95%. Taking B(t — > KinTb> T ) = 0.00273 [22] 
and cr = 87.11 x 10" 6 m [23], we find that ~ 1.4088 and A^ 3 ~ 1.4696. 

4.2 Numerical estimates of the CP asymmetries 

Using the constants for N, N3, and A-D noted above, we integrate Eq. ( [25] ) over 
phase space to obtain dY/dAdx^, dT/dMKn and dT/dM n7T . The results (normalized 
to r tot ) are displayed in Fig. [2j The plots are similar to those in Fig. 2 of Ref. [12] . 
although the agreement is not perfect. We also include a Wess-Zumino contribution 
(not included in the CLEO plots). 
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Table 2: Calculated values for a^} p , a^ Pmod , and Qgp <max - &cP,mod * s computed by 



making the replacement cos if) 



eliminate some of the cancellations that occur upon integration. 



cosijj I in Eqs. (52)-(54). This procedure helps to 

is deter- 



.(*) 



CP, max 



mined by taking the absolute values of the integrands in Eqs. (52)-(54). 



i 


(*) 

a CP 


(*) 

a CP,mod 


(i) 

a CP,max 


CP asymmetry type 


1 


-2.2 x 10~ 5 


-5.2 x 1(T 5 


9.8 x 10~ 4 


Modified rate asymmetry 


2 


7.0 x 10~ 4 


1.0 x 10~ 3 


2.9 x 10~ 3 


Modified rate asymmetry 


3 


2.5 x 10" 4 


6.2 x 10~ 4 


8.3 x 10~ 4 


Triple-product asymmetry 



Having chosen the various coefficients, we can also perform the integrations in 
Eqs. (52)- (54) to obtain the numerical coefficients o!§ P . Numerical values for these 
coefficients are listed in the second column of Table |2j Recall that the actual CP 
asymmetries are obtained by multiplying the a^p by f H lm(r]p) [see Eq. (55) . 

Figure [3] shows plots of the differential asymmetries da^p/dX, with X = Mjc^, 
M„ n , Mktt and cos 6*. In each case, integration over X yields the corresponding 
coefficient a§ P . As is apparent from the figure, each of the asymmetry coefficients 
undergoes considerable cancellation upon integration. Given these cancellations, ex- 
perimentalists may find it advantageous to perform fits to the differential CP asym- 
metries instead of simply measuring the integrated asymmetries. Alternatively, it 
may be possible to achieve larger integrated asymmetries by employing extra weight- 
ing functions in the integration over one or more of the integration variables. As an 
example, we have recomputed the asymmetries with the change cos?/' — ► |cos -01 m 
Eqs. (52)-(54) (as noted above, cosip should be an experimental observable). The 
third column of Table [2] shows the resulting values for the asymmetry coefficients. 
As can be seen, this modification leads to modest increases in the sizes of the coeffi- 
cients. Other weighting functions could also be considered. If a weighting function 
is chosen such that it takes on only the values ±1 over the integration range, the 
largest possible asymmetry coefficients would be obtained by simply taking the ab- 
solute value of the integrand. We have computed the asymmetry coefficients under 
this assumption as well. The results may be found in the fourth column of Ta- 
ble [2] The values in this column represent the maximum values obtainable for the 
magnitudes of the asymmetry coefficients, given the choices we have made for the 
form factors. Comparison of the second and fourth columns in the table illustrates 
the level of cancellation that the integrated asymmetry coefficients have each un- 
dergone. A considerable gain in the magnitude of each asymmetry is possible if an 
appropriate weighting function is adopted. 

A few comments are in order. First of all, we note that the values obtained for 
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Figure 3: Differential asymmetries showing each asymmetry's dependence on the 
respective integration variables. 



the asymmetry coefficients, as well as the shapes of the curves in Fig. [3j depend 
sensitively on the coefficients A-D, N and iV^. We have chosen particular values for 
illustration, but it is assumed that experimentalists would perform more accurate 
measurements of the coefficients A-D in tandem with performing any CP analysis. 
Also, recall that we have assumed that fn is a constant and have thus pulled it 
outside of the various integrations. This may well be a poor approximation, in 
which case the expressions for the differential CP asymmetries would need to be 
modified to include the dependence that fn has on the various variables. Finally, 
we note that more recent analyses use an expression for F s that differs from the 
expression we use [Eq. (59)]. References [TSJ [21] use an expression that is similar 
to Eq. (59), except that it sets a = 1 and A^ 3 = 1, and that it replaces BWk*(Q 2 ) 
and BWk*(s2) by expressions that take into account one or both of the K*'/K*" 
resonances. We have performed an analysis using this modified expression for F^, 
the change affects the asymmetry A CP . Retaining an overall normalization constant 
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and tuning it so that the Wess-Zumino contribution still accounts for approximately 
5% of the experimental branching ratio (A^ 3 ~ 0.4206), we find a^ P ~ 6.8 x 10" 5 
and clqp max ~ 7.1 x 10" 4 . (We do not quote a revised number for a^ PmoA , since 
the replacement cosip — > |cos^| actually makes the magnitude of the asymmetry 
smaller in this case.) Comparison with Table [2] shows that the asymmetries are 
smaller in magnitude in this case. The differential plots are also affected. We do 
not consider results following from this revised expression for F 3 further here, but 
our estimates below could easily be adapted to take this change into account. 

To determine actual CP asymmetries {A^§ P ) from the asymmetry coefficients 
(a,Q P ), we need to know or be able to estimate the quantity /fflm^p) [see Eq. (55) 



Let us begin with a crude estimate by assuming that the NP contribution to the 
width is "hiding" in the experimental uncertainty of the branching ratio. The 
experimental branching ratio determined by BABAR is Bir~ — > K~7t~tt + i> t ) = 
(0.273 ± 0.002 ± 0.009) % |22j; i.e., the experimental measurement has a relative 



uncertainty of approximately 3.4%. A numerical integration of Eq. (49), performed 
under the assumption that only the NP part contributes [i.e., setting Bi=B 2 =B 3 =0 
and _B 4 = B\' = yQ 2 j "hVp / /m r - see Eq. (47)], shows that the experimental un- 



certainty is saturated when \fuVp\ — 17.9. Assuming that rjp is purely imaginary 
yields upper bounds on the magnitudes of the A$ P in the range 3.9 x 10 -4 to 0.012. 
Under the same assumptions regarding fn^p, we also find that the 4p m ax range 
from 0.015 to 0.052. 

The above estimates may be a bit optimistic, although it is difficult to say without 
direct bounds on and r\p. As noted in the Introduction, the CLEO collaboration 
has searched for CP violation in r — > Kkv t \ they have set the following bound on 
the scalar coupling that they denote A [H] , 

- 0.172 < Im(A) < 0.067, at 90% C.L. . (68) 

The coupling A is related to r/s in Eq. ([7]); r/p, however, does not receive a direct 
constraint from this experiment, rjp should scale like m 2 v /m 2 1 due to the Higgs 
propagator (where mw and m# are the W and charged Higgs masses, respectively). 
If the Higgs has electroweak couplings, then it would be reasonable to assume that 
rjp has a magnitude not exceeding unity. At this point we do not have a reliable way 
to estimate One possibility is to infer ju from F 4 using the quark equations of 
motion, although this procedure may well have a large error. As was noted above, 
for this decay has been computed from the perspective of Chiral Perturbation Theory 
in Ref. [IT] . Using the quark equations of motion, one finds \fn\ ~ Q 2 \F±\ /m S) 
leading to an enhancement of fu because of the small strange quark mass. (This 
enhancement would be lost to some degree if the quark mass were replaced by a 
meson mass.) An approximate numerical examination of |i<4| derived from Ref. 



9 We have not updated the expression to account for the possibility of contributions from both 
#1(1270) and ^(1400). 
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shows that it can be of order 1 GeV _1 for some values of Q 2 , s± and S2 (it is also much 
smaller than this for other values of the kinematical variables). A crude estimate of 
the maximum size of \fjj\ would be \fn\ ~ m 2 - x (1 GeV _1 )/m s ~ (1.777 GeV) 2 x 
(1 GeV" 1 )/(0.095 GeV) ~ 30. A more realistic estimate for \fn\ might be in the 
range 1-10. Combining these estimates, we see that |/#Im(7/p)| could be of order 
1-10, leading to a reduction of the possible magnitudes of the CP asymmetries 
compared with our estimates above (for which we assumed |/#Im(7/p)| ~ 17.9). 

5 Discussion and Concluding Remarks 

We have analyzed CP violation in r — > K^tcv t due to NP in the form of a charged 
Higgs boson. Noting that the couplings of a charged Higgs boson to the light quarks 
are suppressed in many models due to the smallness of the light quarks' masses, we 
have observed that CP-odd observables in r — > Kttttv t probe non- "standard" NP 
CP violation. An experimental search for CP violation in r — > KtcttVt would com- 
plement the search for CP violation that has already taken place in r — > Knv T [8]. 
In our notation, r — > Kimv T is sensitive to the coupling r]p, while r — > K-kv t is 
sensitive to r\s- 

We have analyzed four CP-odd observables in r — > K-k-kv t - the rate asymmetry, 
two modified rate asymmetries and a triple-product asymmetry. The rate asymme- 
try is likely to be quite small because it relies on the interference of the SM scalar 
form factor with the NP contribution; thus, we did not make any numerical esti- 
mates for this asymmetry. The modified rate asymmetries and the triple-product 
asymmetry result from the interference of the NP amplitude with the SM contri- 
butions containing the form factors F\ — F$. Adopting a particular model for the 
form factors and making various assumptions, we have estimated the possible sizes 
of the CP asymmetries numerically. In our calculation it was found that each of 
the asymmetries underwent a substantial cancellation upon integration over the 
various phase space variables. Experimentalists may wish to consider differential 
CP asymmetries in order to avoid some of this cancellation. The maximal sizes of 
the three asymmetries (assuming that the cancellations could be avoided by using 
appropriately chosen weighting functions) were found to be in the range 0.015 to 
0.052. These numbers were derived under the assumption that the only constraint 
on the NP contribution is that it is "hidden" in the uncertainty of the branching 
ratio for r — > KirKP T . The maximal magnitudes of the asymmetries decrease if one 
makes more realistic assumptions regarding the hadronic form factor fjj and the NP 
parameter rjp. 

We encourage experimentalists at the B factories to analyze their r data sets 
in the manner that we have described. Future experiments, such as the Super B 
factories, could provide even greater sensitivity to these observables. 

We close with a short comment regarding CP violation in r T — > K^it^K ± v T . In 
principle, this decay mode could be analyzed in a similar manner to what we have 
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described. (See Ref. jl], for example.) One advantage of r — > Ki:Kv T is that there is 
an intermediate pseudoscalar resonance [the 7r'(1300)] that could potentially enhance 
the hadronic current associated with the NP charged scalar exchange. We wish to 
point out what appears to be an error, or an oversimplification, in the literature 
regarding this point. The scalar form factors associated with the n' resonance in 
the r — > 37rz/ r and r — > K-rKv t decays have been written down in Ref. [TU]. The 
expression for the 3tt case seems to be sensible, but the one for the KttK case 
appears to make an unphysical assumption regarding the contributing decay chains. 
In particular, judging from the expression, one of the decay chains would seem to 
have an intermediate p decaying to a K and a 7r. If this is remedied by replacing 
the p by a K*, one finds that none of the decay chains can quite proceed on shell 
(although there is a large uncertainty in the it' mass; furthermore, the ir' does have 
a large width and the decay tt' — > K*K is actually right near threshold). 
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